Astronomy & Astrophysics manuscript no. Fa3 1 1 


February 2, 2008 


(DOI: will be inserted by hand later) 





Far-infrared detection of ^^OH towards Sagittarius B2 
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Abstract. The frequencies and line strengths of the ^113/2 7=5/2-3/2 rotational transition of "OH have been calculated from 
an analysis of its far-infrared laser magnetic resonance spectrum. These results have been used to make the first detection of a 
pure rotational transition of "OH in the ISM. Two resolved components of this transition appear in absorption towards the giant 
molecular cloud Sagittarius B2, which was observed at a spectral resolution of 33 km s"' with the Fabry-Perot mode of the ISO 
Long Wavelength Spectrometer The corresponding transition of "*OH was also observed and its line shape was modelled using 
Hi measurements. The '^O/'^O ratio of 3.5 was then used to compare this with the observed "OH line shape. 
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1. Introduction 

The '^OH molecule was first observed in the ISM towards the 
Galactic Centre by Gardner & Whiteoak J 19761 1 via its mi- 
crowave A-doublet transitions. The strongest of these was then 
observed by Bujarrabal et al. <1983l l towards the giant molecu- 
lar cloud complex, Sagittarius B2 (Sgr B2), with similar obser- 
vations of '''OH and '**OH giving the ratio "^OH/'^OH = 3.6 + 
0.5. Sgr B2 is located ~100 pc from the Galactic Centre and 
emits a bright background ideal for studying absorption fea- 
tures due to the intervening ISM. This spectrum is dominated 
by thermal emission from dust and peaks in the far-infrared 
(FIR) near 80 yum (see Goicoechea & Cernicharo 2001 ). The 
rotational transitions of '^OH lie in this FIR region and have so 
far not been observed in the ISM. 

We report the first detection of a FIR rotational transition 
of '^OH in the ISM. We present accurate frequencies and line 
strengths for the ^03/2 7=5/2-3/2 transition calculated from 
laboratory measurements. We then present observations of this 
transition towards Sgr B2 using the ISO Long Wavelength 
Spectrometer (LWS; Clegg et al. I1996> and model the line 
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shape using the corresponding '^OH line and previously de- 
termined '^0/"0 ratio. 



2. Observations and Data Reduction 

Sgr B2 was observed as part of a wide spectral survey using 
the ISO LWS Fabry-Perot (FP) mode, L03. The whole LWS 
spectral range, from 47 to 197 /vm, was covered with a spectral 
resolution of 30-40 km s"' (see Ceccarelli et al. .2002 for the 
first description of this survey). The lowest energy rotational 
transitions of OH occur in the region 1 19-120 //m. At the res- 
olution of the LWS these transitions have two resolved com- 
ponents due to the A-type parity doubling of each rotational 
level. 

The LWS beam had an effective diameter of 78" at 
119 //m (Gry et al. 2002) and was centred at coordinates a - 
^7h47m2i 75s^ (J = -28°23'14.1" (J2000). This gave the beam 
centre an offset of 2 1 .5" from the main FIR peak, which occurs 
near the radio and mm source, Sgr B2 (M) (Goldsmith et al. 
II992I 1. This pointing was used to exclude the source Sgr B2 (N) 
from the beam. The LWS spectral resolution element at 1 19 jjm 
was «33 km s"'. 

The spectral region 119-120 yum was covered in two sep- 
arate observations using the LWS detector LW2 (ISO TDT 
numbers 50601013 and 50700610) with one further observa- 
tion covering wavelengths below 1 19.74 jum (50800416). Each 
one had a spectral sampling interval of 1/4 resolution ele- 
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ment and 3 repeated scans per point. The three observations 
were processed using the LWS pipeHne version 8 and inter- 
actively calibrated using routines that appeared as part of the 
LWS Interactive Analysis (LIA) package version 10. The dark 
signal, including straylight, was determined as described in 
Polehampton et al. (I^U2t . 

Each observation consisted of a series of FP scans at con- 
secutive grating angles, known as "mini-scans". Due to uncer- 
tainty in the commanded grating angle, the grating response 
profile must be removed from each mini-scan interactively in 
order to recover the true spectral shape (see Gry et al. 2002 (. 
This is particularly important for wide lines that may cross 
one or more mini-scans as it affects the stitching across the 
joins. For weak lines the overlap region can cause extra uncer- 
tainty in the precise line shape if placed unfavourably. A shift 
for each mini-scan was carefully determined using the stan- 
dard LIA tool for interactive processing of FP data, FP_PROC. 
Accurate grating profile shapes were used - these were recov- 
ered from other observations in the Sgr B2 dataset and cor- 
rected for contamination from adjacent FP orders as described 
in Polehampton et al. ( 1200 U . Glitches were removed from each 
observation scan by scan. 

The continuum around the lines was fitted with a polyno- 
mial baseline which was then divided into the data. This effec- 
tively bypassed the large systematic uncertainty in absolute flux 
level (caused by multiplicative calibration steps - see Swinyard 
et al.lI998). The dark signal was accurately determined to be 
a small fraction of the continuum photocurrent and so the fi- 
nal error in relative line depth is dominated by the statistical 
noise in the data. To increase the signal to noise ratio, the three 
observations were corrected to the local standard of rest and 
co-added. The remaining uncertainty in wavelength is less than 
0.004 yum (or 1 1 km s"'), corresponding to the error in absolute 
wavelength calibration (see Gry et al. 2002 1. A more detailed 
description of the calibration and reduction method is given in 
Polehampton ( l2002t . 

3. Laboratory Measurements 

The rotational spectrum of '^OH is more complicated than the 
corresponding spectrum of '^OH and '^OH as there are up to 
6 hyperfine levels associated with the spin of the '^O nucleus 
(7=5/2), each of which is split into two by the 'H hyperfine 
interaction. 

The FIR spectrum of the '^OH radical has been measured 
by laser magnetic resonance (LMR) by Leopold et al. (1987). 
This method involved using a variable magnetic field to tune 
the frequency of the molecular transition into coincidence with 
a fixed frequency laser This meant that zero field frequencies 
were not measured directly in the experiment and were not 
given by Leopold et al. (I1987> . 

In order to predict the relevant frequencies and line 
strengths, we have used a computer program to extrapolate to 
zero magnetic field. This uses the molecular parameters derived 
from the LMR experiment and is described in detail in Brown 
et al. J1978> . This calculation was performed taking account 
of the '^O spin. The spin of causes a further splitting of 
each component into a closely spaced doublet but this was ig- 



Table 1. Calculated hyperfine transition frequencies for 7=5/2- 
3/2 in the ^03/2 state of '^OH. Errors in the last quoted sig- 
nificant figure are shown in parenthesis. An extended version 
of this table for higher transitions up to and including the 
^113/2 7=7/2 level is available electronically at the Centre de 
Donnees Astronomiques de Strasbourg (CDS). 
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nored in the calculation as the doublet separation is well below 
the resolution limit of the LWS FP and does not affect the fi- 
nal line shape. The zero field line frequencies, wavelengths and 
line strengths for each hyperfine component of the 7=5/2-3/2 
transition are shown in Table The line frequencies obtained 
using this technique are accurate to 2 MHz for lines directly 
studied in the LMR experiment. The Einstein coefficient for 
spontaneous emission can be calculated from the line strength 
as described by Brown et al. J1982> . 

4. Results 

Two absorption lines are present in the spectrum of Sgr B2 
which agree well with the calculated wavelengths (Fig. [l) 
and we assign them to the two A-doubling components of 
the 7=5/2-3/2 transition. We are confident that these features 
are real as they occur in all three independent observations. 
Furthermore, the two components agree well in shape as ex- 
pected from all other observed transitions of OH towards this 
source (see Goicoechea & Cernicharo I2002> . They cannot be 
due to "ghost lines" (caused by contamination from adjacent 
FP orders; see Gry et al. i2002i) because there are no strong 
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Fig. 1. Combined data from three L03 observations showing 
the two detected features assigned to '^OH. The velocity scale 
is relative to the strongest components in Table 1: 1 19.6172 fim 
(solid line) and 119.8222 fim (dashed line). The data were 
binned at 1/2 spectral resolution element and the average er- 
ror in each bin across the line is shown. The sharp drop on the 
left hand side of the 1 19.6 /im line is due to the corresponding 
transition in '^OH. 

lines at the wavelengths of neighbouring orders. There is strong 
'^OH absorption next to the detected features but this has a sep- 
aration much less than the offset of nearby FP orders. 

The features show a minimum near the expected velocity 
of SgrB2 (M) (^ +65 km s^'; e.g. Martm-Pintado et al.TWO) 
with additional absorption at negative velocities due to absorb- 
ing clouds corresponding to the main galactic spiral arms cross- 
ing the line of sight (e.g. Greaves & Williams '19941 These 
clouds have a low density and temperature (Greaves .l995j and 
are illuminated by the galactic interstellar radiation field. 

The lines of '*OH and '^OH originating in the ground ro- 
tational state also show this broad structure, although with a 
clearer minimum at the velocity of Sgr B2 (see Goicoechea 
& Cernicharo I2002I I. The two components of the J -5/2-3/2 
transition in '^OH occur at wavelengths of 119.9651 fim and 
120.1718 yum (Morino et al. ll995> and show good agreement 
in shape and depth in the Sgr B2 spectrum. 

5. Line Modelling 

Due to the large errors when directly fitting the weak '^OH 
absorption, the '^OH line was modelled and used with the 
'^O/'^O ratio to predict the '^OH shape. 

Models of photo-dissociation regions show that OH has 
highest concentration at the interface between Hi and H2 
in dense clouds illuminated by UV radiation (Sternberg & 
Dalgarno .l995j . We assume that this is also true in the low 
density clouds towards Sgr B2 and used H i to trace the line of 
sight features. High resolution measurements of the Hi 21 cm 
line have been carried out towards Sgr B2 (M) by Garwood & 
Dickey J1989I I. Ten velocity features are present in their spec- 
trum and we fixed the line widths and velocities to be the same 
for '^OH. The optical depth in each component was adjusted 



Fig. 2. OH absorption line data after co-adding the two A- 
doublet components. The best fit using a model derived from 
H I observations is shown. The data were binned at 1/4 spectral 
resolution element. 

and the spectrum convolved to the resolution of the LWS FP. 
The best fit was found by minimising 

The column density of '^OH in the ground state for each 
component was calculated from the optical depths and lines 
widths. This is a good measure of the total column density 
in the line of sight clouds where only ground state transitions 
have been observed (Goicoechea & Cernicharo i2002.) . The col- 
umn density in the fit summed over the whole line of sight is 
(1.8 + 0.2) X 10'"* cm"^, marginally consistent with the total 
column density found from KAO measurements which gave 
A^(i**OH) > 2 X lO'-* cm-2 (Lugten et al. [T986t . 

In the Sgr B2 envelope, "'OH transitions involving higher 
levels up to 420 K have been observed and radiative trans- 
fer modelling by Goicoechea & Cernicharo (I2002> gave 
A^('**OH) = (6 + 2) X lO'"* cm"2. The column density derived 
from components associated with Sgr B2 in the ''^OH fit was 
(9.8 + 1 .6) X 10'"' cm"^, close to the radiative transfer value. The 
model, results and associated errors and assumptions for '^OH 
have been described in detail by Polehampton (.2002j and wiU 
be presented in a forthcoming paper. 

In order to predict the shape of the '^OH absorption from 
the best fitting model of '**OH, the '^O/'^O ratio was used. 
This has been measured near the Galactic Centre using several 
molecular tracers and these show good consistency (Gardner 
& Whiteoak 1976. Wannier et al. 1976.. Penzias .1981 Guelin 
et al. IT982I Bujarrabal et al. [T983t . The weighted average 
of these values gives '^O/'^O=3.5+0.9, in agreement with 
"*0/'^0 measured throughout the Galactic Disk by Penzias 
(I1981> . However, this is lower than the standard Solar System 
ratio (5.237 + 0.004; Baertschi'T976', Fahey et al.TWt. This 
discrepancy is a long standing problem (e.g. see Prantzos et al. 

As the ratio between "*0 and '^O is not large and both OH 
lines are optically thin, opacity and excitation effects should 
not affect the final ratio. The '^OH column density for each 
line of sight feature was calculated by taking the values from 
the "*OH fit, dividing by 3.5 and splitting up between the hyper- 
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Fig. 3. '^OH absorption line data with model derived from the 
observed '^OH absorption and '^O/'^O ratio of 3.5 (solid line). 
The data were binned at 1/2 spectral resolution element. The 
predicted shape using the Solar System '**0/'^0 ratio of 5.2 is 
also shown for comparison (dashed line). 



fine states within each rotational level (according to the relative 
line strengths of their transitions given in Tabled . This allowed 
the calculation of the optical depth for each hyperfine transi- 
tion and when combined with the line widths and velocities 
from Hi gave a high resolution model of the '^OH absorption. 
The shape has a more complicated structure than for '^OH as 
the hyperfine components have a comparable separation to the 
line of sight features. Figure|3lshows this model (after convolu- 
tion with the FP response profile) with the co-added '^OH data. 
The large number of hyperfine components cause a significant 
modification of the line shape in the convolved spectrum with 
a reduction in the relative absorption of the 65 km s ' feature 
compared with '^OH. Figure|3lalso clearly shows that there is 
too little absorption when the Solar System ratio of 5.2 is used. 

The model fits the '^OH data well in the region of Sgr B2 
but there is extra absorption in the the data near km s"' . This 
may be a real effect as the two resolved A-doublet components 
agree very well in shape. It indicates that the "*0/'^0 ratio 
should be lower than 3.5 in the line of sight clouds, although 
this seems unlikely given the measurements by Penzias fl98H 
that gave a constant ratio throughout the Galactic Disk. It is 
also unlikely that the extra absorption is due to transitions from 
different species as there would have to be two lines with the 
same relative spacing from each '^OH component. 

An instrumental effect that could cause this distortion is 
the mini-scan overlap regions discussed in Sect. 2. Both '^OH 
components occurred in the same position relative to the join 
between mini-scans and so the same effect could have been 
reproduced in each one. Further observations of these lines 
will be necessary to confirm the shape. The accurate wave- 
lengths and line strengths presented here will allow this line 
to be searched for with future instrumentation such as that on 
SOFIA. 
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